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Abstract. The spectral analysis of discretized one-dimensional Schrodinger operators is 
a very difficult problem which has been studied by numerous mathematicians. A natural 
problem at the interface of numerical analysis and operator theory is that of finding finite 
dimensional matrices whose eigenvalues approximate the spectrum of an infinite dimensional 
operator. In this note we observe that the seminal work of Pimsner-Voiculescu on AF embed- 
dings of irrational rotation algebras provides a nice answer to the finite dimensional spectral 
approximation problem for a broad class of operators including the quasiperiodic case of 
the Schrodinger operators mentioned above. Indeed, the theory of continued fractions not 
only provides good matrix models for spectral computations (i.e. the Pimsner-Voiculescu 
construction) but also yields sharp rates of convergence for spectral approximations of op- 
erators in irrational rotation algebras. 



1. Introduction 

In this paper we address the problem of finding numerical approximations to the spectrum 
of a bounded linear operator on a complex, separable Hilbert space. This is a very difficult 
problem in general and there are various ways to attack it as well as various quantities 
that one may wish to use to approximate the spectrum (e.g. pseudospectra or numerical 
ranges). A very natural approach is to start with a given operator T acting on a Hilbert 
space H and "compress" it to a finite dimensional subspace of H. This compression is a 
finite dimensional matrix whose eigenvalues can hopefully be numerically computed and one 
further hopes that these eigenvalues will somehow approximate the spectrum of the infinite 
dimensional operator T. This approach has been studied by numerous authors (see, for 
example, [6], [2] and the references therein) and a number of interesting results have been 
obtained for large classes of operators. 

Though the present work grew naturally out of the author's own study of the method 
described above (cf. [4, Section 6]) it turns out that an absurdly abstract approach to nu- 
merical approximations of spectra actually proves useful in at least one important case. The 
main result of this paper simply observes how some very specialized and technical work in 
C*-algebra theory can be used to get excellent approximations to spectra of some important 
operators. In particular, this applies to discretized one-dimensional Schrodinger operators 
with quasiperiodic potential. See [1] for a nice treatment of the spectral theory of Almost 
Mathieu operators (up to 2001) and [3] (and its references) for a recent contribution to the 
more general case of quasiperiodic Schrodinger operators. 

When presented abstractly the basic idea of this note becomes quite simple, but will require 
a bit of C*-algebra theory to explain. The first basic fact we need is that all finite dimensional 
C*-algebras have a special form - they are just finite direct sums of finite dimensional matrix 
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algebras. We will also need the definition of an AF ("approximately finite dimensional") 
algebra; A C*-algebra C is AF if there exist finite dimensional C*-subalgcbras Ci C C2 C 
■ ■ ■ C C whose union is dense in C. The observation below is the trivial, but key, idea of 
this paper. 

Observation: If T e B{H) is an element in an AF C*-subalgebra of B{H) then there 
exist finite dimensional matrices whose spectral quantities (e.g. pseudospectra or numerical 

ranges) approximate those of T. 

Indeed, if T e C C B{H) with C an AF algebra then we can find operators T„ such 
that yr — T„|| — > (hence spectral quantities of T„ are close to those of T) and each T„ is 
contained in a finite dimensional C*-algebra (hence can be identified with a finite dimensional 
matrix). There is one subtlety that must be mentioned. Namely, the actual spectrum of T„ 
need not be close to the actual spectrum of T (unless T happens to be normal) but singular 
values, numerical ranges and pseudospectra will be close. Moreover, the actual spectrum 
of T is equal to the intersection of all its pseudospectra and hence we still get reasonable 
approximations to the spectrum of T by looking at pseudospectra of the T„'s. In the case 
both T and T„ are normal (e.g. self-adjoint Schrodinger operators) then the actual spectra 
are close in a very strong sense; the Hausdorff distance between their spectra is bounded by 
IIT - T II 

II 

AF-embeddability (i.e. studying which operators belong to AF algebras) has been consid- 
ered by numerous authors but most of this work is too general to be of much help in finding 
explicit matrices which approximate given operators. For example, it is an easy consequence 
of the spectral theorem that every normal operator is contained in an AF-subalgebra of 
B{H) but this gives no clue how to actually find the right matrix approximations. However, 
returning to the seminal paper on AF-embeddability (cf. [8]) one finds that not only arc ex- 
plicit matrix models provided (for operators in irrational rotation algebras) but computable 
and sharp rates of convergence can be proved. 

In section 2 we will review the necessary aspects of the theory of continued fractions. We 
also state a few simple lemmas that will be needed later. 

In section 3 we review the construction of Pimsner-Voiculescu which provides explicit 
matrix models for AF-cmbcddings of the irrational rotation algebras. We will not reproduce 
the proofs of any estimates as these can be found either in the original paper [8] or in [5, 
Chapter VI]. 

In section 4 we state and prove the main result for spectral approximations of operators 
in irrational rotation algebras. The proof really amounts to working out a few estimates as 
the hard part was accomplished more than 20 years ago in [8]. 

Finally, in section 5 we show how results of Haagerup-R0rdam can be used to give general 
"one-sided" rates of convergence. Though the results of this section are not as good as 
previous sections they have the advantage of always being numerically implementable. 



2. Preliminaries 

In this section we will review some basic facts about continued fractions and point out a 
few technical (but very elementary) inequalities that we will need later. Though there are 
a number of excellent books on continued fractions, a classic reference is [7]. The facts we 
will need, however, can be found in any book on the subject. We also state the definition of 
pseudospectrum at the very end of this section. 
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Given an irrational number 9 e (0, 1) there exist unique positive integers ai, 02, . . . such 
that ^ 

9 — hm [ai , . . . , a„] = hm 



a2- 



Note that if 9 is given then one can compute the integers a„ via the following recursive 
procedure: Define 9i = ^, ai — [9i\, 9n+i — g , dn+i — L^n+iJ, where [x\ is the integer 
part of a real number x. As is customary, we use p„, Qn to denote the numerator and 
denominator, respectively, of the n*'^ convergent 

r 1 _ 1 -P^ 
ai, . . . , a„ — ■ 1 — — . 

«1 + a, + ^— 'i- 

'■■+^ 

an 

The PnS and q^s satisfy the recursive relations 

Pn = (^nPn-l + Pn-2) In = O.nQn-l + Qn-2, 

where po = 0, pi = 1, qo = 1 and qi = ai. The following remarkable fact can be found in 
any book on continued fractions. 

Theorem 2.1. For every n one has 1^ — — I < — ^ — < i—Y. 

We now record a few trivial lemmas that will be used in later estimates. 

Lemma 2.2. IfF{k) is the k^^ Fibonacci number (where F(0) = 1 = F{2) = 2, etc.) 

then for each pair of positive integers n and k we have qn+k > F{k)qn. 

Proof. Use induction and the recursion formula. □ 

00 

Lemma 2.3. E ^ < 

k=Q 

Proof. Induction shows that F{k) > and hence the result follows by taking recip- 

rocals and applying the formula for a geometric series. □ 



oo 

Lemma 2.4. E — < -^r^- 



fe=0 

Proof. Apply the two previous lemmas. □ 

We close this section with the definition of pseudospectra. 
Definition 2.5. For e > 0, the e-pseudospectrum of T e B(H) is 

(t(^)(t) = {A e c : ||(A-r)-i|| > -}, 

e 

where = cxd if X G B{H) is a non- invert ible operator. Note that the usual spectrum 

o^{T) is contained in a^^\T) for every e > and, moreover, 

£>0 

See http:/ /web. comlab.ox.ac.uk/projects/pseudospectra for more, including software 
for computing pseudospectra. 
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3. The Pimsner-Voiculescu Construction 

In this section we state the technical aspects of the Pimsner-Voiculescu construction that 
we will need. We see no reason to reproduce proofs as all the details can be found in [5, 

Sections VI. 4 and VI. 5]. 

Throughout this note, Af^.(C) will denote the k x k complex matrices. Let < p < q he 
integers, cu = e '^^'^ and define matrices up,vp G Mg(C) as follows. 
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If we are given an irrational number 9 G (0, 1) and we let ^ denote the n*^ convergent (as 
in the previous section) then to ease notation we will let 6n = —, uq^ ~ uvn and vq^ — Veil. 

Theorem 3.1. (Pimsner-Voiculescu) For each natural number n there is a *-monomorphism 

■■ Mg^_^{C) © M,„(C) ^ M,„(C) © M,„^,(C) such that 

II / N II 27r 

and 

\\<Pn{ver,-, ® vej - ve^ © ve^+, \\ < + — . 

Qn-l Qn 



4. Computing Spectra in Irrational Rotation Algebras 

For an arbitrary number 9 e (0, 1) there is a universal C*-algebra, denoted by Ag, which 
is generated by two unitaries Ug, Vg subject to the commutation relation 

Universality means that iiUg, Vg are any other unitary operators satisfying the same relation 
then there exists a surjective *-homomorphism Aq — > C*{Ue, Vq) such that Uq i— > Uq and 
Vq ^ Vg. One crucial fact which we will need is that if 9 is irrational then Ag is a simple 
C*-algebra (i.e. has no non-trivial, closed, two-sided ideals) and hence any *-homomorphism 
will be injective in this case (cf. [1, Theorem 1.10]). 
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In the theorem below dni-,-) denotes the Hausdorff distance between two compact subsets 
of the complex plane: 

dni^, A) = max{sup d{a, A), sup d{X, E)}, 
o-eE AeA 

where d{a, A) = inf^gA \(J — A|. 

Theorem 4.1. Let 9 G (0, 1) he irrational, 9n = tt±ii/9±i & C be four complex numbers 
and M = max{|Q:±i|, // we let Hg = aiUg + a-iUg + (3iVe + and h^^ — 

oiiue^ + a-iul^ + I3ive„ + then: 

(1) (cf. Definition 2.5) For each e > we have 



(e+2£„) 



{her.: 



where = 204M(^ + 
(2) (Normal Case) If ho^_j^,hg^ and Hg all happen to be normal operators (e.g. self- 
adjoints) then 

dnHhe^.,) U (7(/i,J, <T{He)) < 204M(^ + -). 

Qn-l Qn 

Proof. It follows from the Pimsner-Voiculescu construction that there is an AF algebra 21^ 
and a sequence of *-monomorphisms : My^_^(C) © Mq^(C) ^ Sl^ with the property that 

27r 

Qn-l 

and 

71" 47r 

||V'n(^^0„_i © vej - i^n+i{ver. ® ve^+i)\\ < + — • 

Qn-l Qn 

Due to the exponential growth of the QnS we have that both {i/j^ ©we„ ) } and {tpn (^^n-i © 
vo^)} are Cauchy sequences of unitaries and hence converge to some unitaries in Stg which 
satisfy the same commutation relation as Ug and Vg. By universality there exists a (nec- 
essarily faithful - by simplicity) *-homomorphism ^ : Slg with the property that 
i'niuen-i © %„) "^{Ue) and ilJn{ve„-i ® vgj "^{Ve). More importantly, we have the 
following estimates: 



fc=0 

27r 



fe=o ^'^+'^-1 



= 27r( + — +>" ) 

1n-l Qn ^ Qn+k+1 

^ , 1 1 , 47rV5 1 

< 27r( + — ) ■ 



Note that we applied Lemma 2.4 to get the last inequality. A similar argument shows that 
\\-i'{Ve)-Mv0n-^®veJ\\<-— + 



Qn-l Qn a/5 — 1 Qn+1 
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With these estimates in hand it is immediate that 

\\^iHo)-Uher.-.®heJ\\ < {\a,\ + \a,\){2ni^ + -) + '^'''^ ^ 



' Qn-l Qn — 1 Qn+1 

/I ^ I I /I IX / ^ Stt IOttV^ 1 X 

+ m + m){ — + — + ^^-r — • 

Qn-l Qn V5-l9n+l 

The estimate above, though a bit ugly, is the one which should be used in practice as the 
constants are smaller than those claimed in the theorem. To get the cleaner statement of 
the theorem we first note that the recursion formula Qn+i = + implies 

111 

< + -, 

Qn+l Qn-l Qn 

and thus 

\\^{He)-Mher.-.®heJ\\ < 2m( 2n + + -) 



V^S - 1/ Qn-l Qn' 

^^(^ IOttVsX, 1 1, 
+ 2M 57r+ ^ + — 

V V5 - 1 / Qn-\ Qn 

= 2M 77r+ + — 

V V5 - 1/ Qn-\ Qn 

,^3^5-1, 1 1, 
= 147rM— ^ ( \ ) 

V5 - 1 Qn-\ Qn 

< 204M(^ + — ). 

Qn-l Qn 

Now we quote some general results concerning approximation of spectral quantities. Given 
two operators S,T e B{H) one has: 

(1) (of. [6, Theorem 3.27]) For every £ > 0, a^'\S) C a(^+ll^-^ll)(T) C a^'+-'\\^-n)[S). 

(2) If both 5" and T happen to be normal then dH{(T{S),a{T)) < ||5 - T||. 

With these general facts and the norm estimates above one easily completes the proof by 
observing that the spectral quantities of © T (i.e. pseudo or actual spectrum) are just the 
union of the spectral quantities for S and T. □ 

Remark 4.2. Almost tri-diagonality. Note that h^^ is "almost" a tri-diagonal matrix. It is 
non-zero in the upper-right and lower-left entries but otherwise is tri-diagonal. In particular, 
these matrix models are sparse. 

Remark 4.3. Discretized Schrddinger operators with polynomial potential. 

The theorem above is stated for a class of operators that includes all of the examples dis- 
cussed in [1] (i.e. Almost Mathieu operators and a few non-self- adjoint examples). However, 
similar results hold for a broader class of examples that includes all the one-dimensional 
discretized Schrodinger operators with polynomial potential. In fact, if R{X, Y) is any poly- 
nomial in non-commuting variables X and Y then 

W^iRiUe, Ve)) - ^„(i?(«.„_„ t;.„_J © R{ue^,veM = + -). 

Qn-l Qn 

Unfortunately, the constants involved get quite complicated quite quickly but they can be 
explicitly worked out for a given polynomial should one so desire. Having this norm estimate 
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one gets similar convergence results for spectral quantities in this setting. Finally we should 
point out that the spectral quantities of any operator in an irrational rotation algebra (e.g. 
Schrodinger operators with arbitrary quasiperiodic potential) can, in principle, be approxi- 
mated with this method - first approximate by a polynomial and then proceed - however, 
it seems impossible to control rates of convergence in this generality. 

Remark 4.4. These rates of convergence are sharp. 

To see this we consider the case that Hq — Uq and 9 is an irrational number for which the 
a„'s are uniformly bounded. Then dH{a{H0),a{h0^_^) U a{h0j) < 292(^ + ^). It IS easy 

to check that — — h — < and hence if one could prove that dH(o'(Ha),a(he .) U 

q„-l qn — q„-i+qn ^ n\ \ v/, \ v„-i/ 

(^(hej) = + ^) then it would follow that dH{(T{Hg),a{h0„_J U a{h0j) = o( g^_^+gJ - 

However this is impossible as it would imply that the Lebesgue measure of (T{Hg) is zero 
(which it isn't: cr{U0) = (A G C : |A| = 1}) since there are at most g^-i + Qn points in 

Remark 4.5. Implementability. It is, of course, true that for some irrational numbers the 
strategy proposed in this paper for approximating spectra can't reasonably be implemented. 
For example, it can happen that qi is very small while §2 is far larger than any computer 
can handle (just take a2 to be huge). On the other hand, the set of irrationals for which the 
associated sequences {a„} stay bounded is a dense set in (0, 1). For example, all irrationals 
which are the roots of quadratic equations have this property; hence i/r for any rational 
number and these are already dense in (0, 1). In other words, for a dense set of irrationals this 
program is reasonable and hence the question becomes how spectra behave as the parameter 
6 changes. But it is well known that spectra of operators in irrational rotation algebras vary 
continuously in 9 (see the next section) and hence the program suggested in this paper is 
about as practical as could be hoped for. 

Remark 4.6. Uniqueness of convergents. Another basic, but remarkable, property of con- 
tinued fractions is that if 9 is irrational and \9 — -\ < then ^ is necessarily one of the 
convergents in the continued fraction decomposition of 9 (i.e. if ^ = [oi, 02, . . .] then there is 
an n such that | = [oi, . . . , a„] = - see [7, Theorem 184]). The point is that if one starts 
with a rational | and only looks at the irrationals 9 which are a distance at most ^ away 
from - then the spectral quantities of hp will provide good approximations to those of Hg 

Q q 

(since | necessarily arises in the continued fraction decomposition of all such 9 and hence 
the theorem above applies). Of course, rates of convergence become trickier in this setting 
but our only point is that one can simultaneously approximate a small interval of operators 
with a single matrix model. 

5. General One-sided Rates of Convergence 

In our final section we will record a consequence of the following theorem of Haagerup and 
R0rdam (cf. [1, Corollary 3.5]). Among other things this shows that if one fixes a polynomial 
in the canonical generators then the spectral quantities vary continuously in 9. 

Theorem 5.1. For every 9,9' e [0,1) there exists a Hilbert space H and injective *- 
homomorphisms irg : Ag ^ B{H), irgi : Agr — > B{H) such that 

he'iUg,) - 7^6{Ug)l \\T:g\Vg,) - 7r,(Ve)|| < ^^Qti\9-9'\. 
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As in the last section we will restrict attention to operators of the form Hg — aiUo + 
o:_iUg + + P-iVg as the estimates are cleaner in this setting. However, the same 
techniques handle more general polynomials as well. In the theorem below we only get "one- 
sided" rates of convergence (as opposed to estimates on the Hausdorff distance). The modest 
advantage, however, is that this strategy is implementable for any irrational number. 

Wc find the following notation convenient: For compacts A, E C C and S > we write 
A S if for each A G A there exists cr e E such that \X- (t\ < 5. Note that dniK S) < 5 
if and only if A S and S A. 

Theorem 5.2. Let 9 G (0, 1) be irrational and for each n G N choose an integer pn G 
{0, 1, . . . ,n — 1} such that \0 — ^\ < 2n' operator polynomial as above, hn = 

aiUEn +a-iu*pn + PiVpn + (3-iv*pn is the associated matrix model and M — max{|a-|-i|, |/3±i|} 

n n ^ 

then the following statements hold: 

(1) For every £ > 0, 

where Ci < 36MV^. 

(2) If Hg happens to be normal then 

a{hn) a{He), 

where Ci < 36MV37r. 

Proof. To prove this result we first observe that all of the spectral quantities associated to /i„ 
are contained in the corresponding quantities for Hpn — aiUpn + cc-if/l^ + PiVsn + P^iV^ 

n 71 n 

(since universality gives us a *-homomorphism taking Hpn 1— > hn)- The Lip^/^-continuity 
theorem of Haagerup-R0rdam then gives upper bounds on the distance between the spectral 
quantities of Hq and Hpn which completes the proof. □ 

n 

Remark 5.3. Hausdorff Convergence. It can be shown that dH{(J^^\hn),a^^\He)) (for 
every e > 0). The proof is similar to the proof of [4, Theorem 3.5] and hence will be omitted. 
However, it does not seem possible to control the rate of convergence in the Hausdorff metric. 
For example, it is certainly not the case that rates of convergence are controlled by \0 — ^{'^ 
for any a > because many irrational numbers (e.g. transcendental numbers) have extremely 
good rational approximations. To be more precise, it can happen that a particular 6 has 
the property that there are infinitely many co-prime solutions to the inequality 1^ ~ ^1 < ^ 
where k is some fixed natural number. If the Hausdorff distance between spectra could be 
bounded above by 16* — ll'^ then we would also get an upper bound of (^)'"^. However, when 
ka > 1 it would then follow that the spectrum of Hg has Lebesgue measure zero since there 
are at most q eigenvalues for the matrix hp . Since these operators need not have spectra of 

Q 

measure zero it follows that an upper bound of the form \9 — ^{"^ is impossible. 
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